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Abstract 

Bid-ask spread is taken as an important measure of the financial market liquidity. In this 
article, we study the dynamics of the spread return and the spread volatility of four liq- 
uid stocks in the Chinese stock market, including the memory effect and the multifractal 
nature. By investigating the autocorrelation function and the Detrended Fluctuation Anal- 
ysis (DFA), we find that the spread return is lack of long-range memory, while the spread 
volatility is long-range time correlated. Moreover, by applying the Multifractal Detrended 
Fluctuation Analysis (MF-DFA), the spread return is observed to possess a strong multi- 
fractality, which is similar to the dynamics of a variety of financial quantities. Differently 
from the spread return, the spread volatility exhibits a weak multifractal nature. 
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1 Introduction 



Financial dynamics has attracted a central interest of physicists in recent years yj- 
ls|]. The financial market has accumulated a large amount of historical data in the 
past years, which provides an abundant data resource for the empirical study. Some 
scaling behaviors have been revealed to be universal for different markets [[ihQ], 
e.g., the well known "inverse cubic law" of the probability distribution of the price 
return. The so-called "volatility clustering" is also found for a variety of financial 
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markets [0, [kJ, which refers to the long-range time correlation of the absolute 
value of the price return. Differently from those universal behaviors, unique char- 
acteristics are found for several emerging markets. A typical example is the "anti- 
leverage effect" of the Chinese stock market 111 in . which indicates an anti-effect 
from the "leverage effect" of most mature financial markets in the return-volatility 
correlation lll2n . On the other hand, the financial markets can be taken as a com- 
plex system composed of many elements, with the scaling behavior emerging from 
the many-body interactions. Different models and theoretical approaches have been 
developed to describe financial markets I113l-ll7ll. 



Within this framework, the so-called bid-ask spread is served as an important indi- 
cator to quantify the financial market liquidity and efficiency. Most modern finan- 
cial markets are order-driven markets, which adopt the continuous double auction 
mechanism. There are two basic types of orders, the market order and the limit or- 
der. For the market order, trade is carried out immediately when receiving the order. 
For the limit order, order is stored, and it waits for a best price for trading. The buy 
limit order is called bid, and the sell limit order is called ask. The bid-ask spread 
refers to the price difference between the lowest ask price and the highest bid price. 
Generally speaking, the lower the spread, the smaller the transaction cost, and also 
a higher liquidity of the market. 



So far, statistical properties of the bid-ask spread of the limit-orders have been 
widely studied lll8l-l2lll. Empirical study shows that, the probability distribution of 



the bid- ask spread obeys a power law behavior, with the exponent around 3.0 119, 



20I1 . Long-range time correlation of the bid- as k sp read is also revealed for different 



markets, including the Chinese stock market |U21]. In addition, differently from 
the multifractality of a variety of financial quantities, the bid-ask spread is reported 
to be mono-fractal for the Chinese stock market 112 lh . On the other hand, a lot 
of work contributes to the order flow dynamics, and the microscopic mechanism of 
price formation Ill8l- l20l l22l-l27n. However, the dynamics of the spread return and the 
spread volatility has not been analyzed or reported in detail, to our knowledge. The 
spread return is defined as the difference of bid-ask spread between two successive 
time point, and the spread volatility is defined as the magnitude of the spread return. 
In fact, in finance, one concerns more on the price return or the price volatility, 
rather than the price itself. Similarly, it is essential to investigate the bid-ask spread 
return and the spread volatility dynamics. 



With this incentive, in this paper, we study the dynamics of the spread return and the 
spread volatility. Based on the tick-by-tick data of four liquid stocks in the Chinese 
stock market, we investigate the memory effect of the spread return and the spread 
volatility by applying the autocorrelation function and the Detrended Fluctuation 
Analysis (DFA). Multifractal nature of the spread return and the spread volatility is 
also revealed, by employing the Multifractal Detrended Fluctuation Analysis (MF- 
DFA). 
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The remainder of this paper is organized as follows. In the next section, we present 
the datasets we analyzed. In section 3, we introduce the definition of the bid-ask 
spread, the spread return, and the spread volatility, and then we focus on the mem- 
ory effect of the spread return and the spread volatility by applying the autocorre- 
lation function and the DFA method in section 4. In section 5, we emphasize the 
multifractal nature of the spread return and the spread volatility. Finally comes the 
conclusion. 



2 Datasets 



The Chinese stock market is an order-driven market based on the continuous double 
auction 112 ll |28L 12911 . Before July 1, 2006, only limit orders were permitted to sub- 
mit in the order placement. The Chinese stock market is composed of two stock ex- 
changes, the Shanghai Stock Exchange (SHSE) and the Shenzhen Stock Exchange 
(SZSE). The SHSE was established on November 26, 1990, and the SZSE was 
established on December 1, 1990. The dataset we analyzed is based on the tick-by- 
tick data of 4 liquid stocks listed on the SHSE and the SZSE, i.e., the Shenzhen 
Development Bank Co., Ltd (SDB) stock, the Guangdong Electric Power Devel- 
opment Co., Ltd. (GEP) stock, the Datang Telecom Co., Ltd (DTT) stock, and the 
Shanghai Lujiazui Finance and Trade Zone Co., Ltd (SLFT) stock. The time scale 
of the datasets covers 3 whole year from 2004 to 2006, and the time resolution of 
the data record is about 15 seconds on average. 



3 Definition of the Spread Return and the Spread Volatility 



Different definition of the bid-ask spread has been proposed II 1 8l |20L l24l |30h35I1 . In 
our study, we define the bid-ask spread S as the price difference between the lowest 
ask price and the highest bid price, which reads, 



S(t') = a(t')-b(t'), 



(1) 



where a(t') and b(t') arc the lowest ask price and the highest bid price at time t', 
respectively. To quantify the spread S at a same time scale, we define a rescaled 
spread S in a time interval At as 1180, 



JV 



iV i=l 



(2) 



where N is the total number of transactions in the time interval At. In this article, 
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we study the time scale At = 1 minute. The original spread return is defined as the 
spread change between two consecutive values of the rescaled spread, 



R(t') = ln( 



m 
sit- 



es) 



The original spread volatility is then defined as the absolute value of the original 
spread return, i.e., V(t') = \R{t')\. 



It has been reported that the bid-ask spread shows similar characteristics as a num- 
ber of financial quantities, such as the long-range time-correlation of the spread 
series. However, the spread is also found to present some unique feature. For ex- 
ample, the spread exhibits a mono-fractal characteristic 112 ill , which is different 
from the multifractality of most financial quantities. What feature could the spread 
return and the spread volatility present? In the following section, we investigate the 
memory effect of the spread return and the spread volatility by the autocorrelation 
function and the DFA method. 
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Fig. 1. The autocorrelation function Aj^(t) of the original spread return R(t') is displayed. 



4 Memory Effect 



4. 1 Autocorrelation 



It is well known that the price return exhibits a short-range time-correlation, while 
the price volatility is long-range correlated. To detect the memory effect of the 
spread return and the spread volatility, we investigate the autocorrelation function, 
which is defined as, 
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Az(t) = [(R(t')R(t + t'))-{R(t'm/a 



(4) 



Ay(t) = [(V(t')V(t + 0> - (VW) 2 ]/4, 



(5) 



where a\ = (R(t') 2 ) - (R(t')) 2 , and o 2 ~ 
an average over t' . 



{V{t') 2 ) - (V(t'))\ and the (■ ■ ■} takes 
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Fig. 2. The autocorrelation function A~(t) of the original spread volatility V(t') is dis- 
played on a log-log scale. 

In Fig. 1, the autocorrelation function of the original spread return R(t') is dis- 
played for the SDB, the GEP, the DTT, and the SLFT stock, respectively. It is 
observed that, for the four stocks, the autocorrelation of the spread return fluctuates 
around 0, which indicates a short-range time-correlation. 

Fig. 2 shows the autocorrelation function of the original spread volatility V(t'), 
where a periodic structure is observed. Such kind of periodic structure refers to 
the so-called intraday pattern, which has been found in a variety of high-frequent 
financial series D6H41I1 . including the bid-ask spread 112 ll |42[ 14311 . In the Chinese 
stock market, there is an active opening call auction time period from 9:15 to 9:25 
a.m., when the buy and sell orders are aggregated to match. The intraday pattern 
could be dominated by the active call auction at the opening time. The intraday 
pattern strongly influences the financial dynamics, therefore it requires an efficient 
method to remove the intraday pattern. A $-time scale method has been applied 
to remove the intraday-pattern of a FX market IBal . In our study, we follow the 
method introduced by Ref . 



In the Chinese stock market, there are 240 consecutive working minutes in a trading 
day. We then segment the dataset of each trading day into 240 consecutive 1-min 
intervals. We define the intraday pattern of the spread return and the spread 
volatility My at time t' in these 240 working minutes as, 
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1 N d 



M R(t') = w i: R i(t% (6) 



1 N d 



My(t') = Wd Y,V^')i (7) 



1=1 



where i is the i th trading day, Nd is the number of trading days. The adjusted spread 
return R and the adjusted spread volatility V is then defined as, 



R(t') 



M~ R (t>) 



(8) 



V(t') 



V(t') 

My(t>) 



(9) 



The autocorrelation function A v (t) of the adjusted spread volatility is defined in a 
similar way as that of the original spread volatility, by replacing the original spread 
volatility V with the adjusted volatility V. After removing the intraday pattern, the 
autocorrelation function of the adjusted spread volatility is shown in Fig. 3, where 
most of the periodic peaks are eliminated. The Ay(t) is found to decay slowly, 
which suggests a long-range time correlation of the spread volatility. However, the 
autocorrelation dynamics shows a large fluctuation, therefore it is difficult to detect 
the decay exponent of the A v (t). 




t (min) 
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Fig. 3. The autocorrelation function Ay(t) of the adjusted spread volatility V(t') is dis- 
played on a log-log scale. 
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4.2 DFA Method 



For the non-stationary time series, a widely adopted measurement of the time- 



correlation is the DFA method 11441 14511 . which has been successfully applied to 
detect the time correlations in various physical systems. We then introduce the DFA 
method. 

Considering a fluctuating time series A(t'), one can construct a new time series by, 

f 

B(t')= ~ A ave ], (10) 

t"=l 



where A ave is the average of the A(t') in the total time interval [1, T]. By dividing 
the total time interval into windows N t with a size of t, and linearly fit B{t') to a 
linear function B t {t') in each window, the DFA function of the k th window box is 
defined by, 



fkitf = - £ [Bit') - B t {t')]\ (11) 
1 f=(fc-i)t+i 

The overall detrended fluctuation is then estimated as, 



i N t 

^) 2 =atEM' (12) 
iv < k=l 

In general, F 2 (t) scales by a power-law F 2 (t) ~ t H . If 0.5 < H < 1.0, A(t') is 
considered to be long-range correlated in time; if < H < 0.5, A[t') is temporally 
anti-correlated; H = 0.5 corresponds to the Gaussian white noise, while H = 1.0 
suggests the 1/ / noise of the A(t') series. If H is bigger than 1.0, the A(t') series 
is taken as an unstable series. 

We apply the DFA to the adjusted spread return and the adjusted spread volatility. In 
Fig. 4, the DFA function of the adjusted spread return is shown on a log-log scale. A 
power-law increase is observed, with the exponent is measured to 0.49 for the SDB, 
0.48 for the GEP, 0.49 for the DTT, and 0.51 for the SLFT. The DFA exponents of 
the four stocks are all around 0.50, which indicates a Gaussian white noise of the 
adjusted spread return, i.e., no memory in the spread return series. The results of the 
DFA method are consistent with those obtained from the autocorrelation functions 
of the spread return. The absence of memory in the spread return series is similar as 
that of the price return. Differently from the spread return, the DFA of the adjusted 
spread volatility presents a cross-over phenomenon for all the four stocks, as shown 
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Fig. 4. The DFA function F(t) of the adjusted spread return R is displayed on a log-log 
scale. 



in Fig. 5. For the forepart time scale, the DFA exponents are measured to be 0.68 
for the SDB, 0.73 for the GEP, 0.65 for the DTT, and 0.66 for the SLFT. The value 
of the four exponents is ranged inbetween (0.5, 1), which indicates a long-range 
correlation of the spread volatility. For the afterward time scale, DFA exponents 
are measured to be 1.00 for the SDB, 1.07 for the GEP, 1.11 for the DTT, and 
0.81 for the SLFT. Except for the stock SLFT, the other 3 stocks present a 1/ f 
noise characteristic. Such kind of cross-over behavior of the spread volatility is 
also similar to the price volatility dynamics. 

To this end, the spread return is found to be lack of long-range memory, whereas 
the spread volatility exhibits a long-range time correlation. The memory effect of 
the spread return and the spread volatility is quite similar to that of the price return 
and the price volatility, respectively. 
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Fig. 5. The DFA function F(t) of the adjusted spread volatility V is displayed on a log-log 
scale. 
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5 Multifractal Analysis 



The multifractal analysis is widely used to evaluate the risk of the financial market 
Il46ll . also to quantify the market inefficiency I147n . Multifractal characteristic has 
been found in many quantities of financial markets, such as the price volatility, and 
the trading volume Il48l-l50ll. However, differently from most of the financial quan- 
tities, the bid-ask spread has been reported to be mono-fractal. Could the spread 
return and the spread volatility possess the multifractal nature? In our study, we 
apply the MF-DFA to reveal the multifractal feature of the spread return and the 
spread volatility. The MF-DFA is a generalization of the original DFA method. It 
considers different orders of the detrended fluctuations, which reads, 



1 N t 

m = { 7 rY.imn 1/q , 03) 
iv * k=i 

where q is the q t h order of the detrended fluctuation. If q — 2, the MF-DFA then 
reduces to the original DFA method. The q can take any real number except q = 0, 
where F (t) reads by, 



i Nt 

F (t) = exp{- (14) 

* k=l 

Generally speaking, the MF-DFA function F q {t) scales with the window size t by, 



F q {t)~t h ^\ (15) 

where h(q) is the MF-DFA exponent. For a large absolute value of q, the F q (t) 
shows a large fluctuation with t, due to the finite size of the investigated financial 
system. Therefore, in our study, we investigate the multifractality of the spread 
return and the spread volatility with q E [—6, 6]. 

In Fig. 6, the F q {t) of the adjusted spread return R and the adjusted spread volatility 
V is shown for the SDB stock and the DTT stock. It is observed that, the F q (t) of 
both the adjusted spread return R and the adjusted spread volatility V shows a 
power law scaling for the investigated value of q, except for a two- stage power law 
scaling in the F q {t) of the adjusted spread return R of the DTT stock at q = —2, —4. 
The cross-over behavior of the DTT stock might be from fluctuations induced by 
the big value of q. From Fig. 6 (a), the exponent of the F q {t) is estimated to be 
1.20, 1.10, 0.80, 0.49, and 0.35 for the q = -4, -2, 0, 2, and 4 for the SDB stock, 
respectively. That is to say, the increasing trend of the F q (t) with t becomes more 
slowly as q increases. Similar behavior is also found for the DTT stock. It indicates 
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Fig. 6. (a) is the MF-DFA function F q (t) of the adjusted spread return for the SDB stock, 
with the exponents measured to be 1.20, 1.10, 0.80, 0.49, and 0.35 for q = -4, -2, 0, 2, 
and 4, respectively, (b) is the F q {t) of the adjusted spread volatility for the SDB stock, 
with the exponents measured to be 0.77, 0.75, 0.74, 0.73, and 0.71 for q = -4, -2, 0, 2, 
and 4, respectively, (c) is the MF-DFA function F q (t) of the adjusted spread return for the 
DTT stock, with the exponents measured to be 0.59, 0.50, and 0.44 for q = 0,2, and 4, 
respectively. For q = —4 and —2, it shows a cross-over behavior, with exponent measured 
to be 1.21 and 0.96 for the forepart time scale, and 0.74 and 0.65 for the afterward time 
scale, respectively, (d) is the F q (t) of the adjusted spread volatility for the DTT stock, with 
the exponents measured to be 0.66, 0.66, 0.65, 0.64, and 0.63 for q = -4, -2, 0, 2, and 4, 
respectively. 



a multifractal nature of the spread return. For the adjusted spread volatility, the 
exponent of the F q (t) is estimated to be 0.77, 0.75, 0.74, 0.73, and 0.71 for q = -A, 
-2, 0, 2, and 4 for the SDB stock, respectively. The exponent of the F q {t) decreases 
very slowly with q. Also, similar behavior is observed for the DTT stock. It suggests 
a weak multifractal feature of the spread volatility. 

To further understand the multifractal nature of the spread return and the spread 
volatility, we investigate the MF-DFA exponent h(q) on q, the scaling exponent 
function r(q) on q, and the multifractal singularity spectrum f(a) on a, for the 
adjusted spread return and the adjusted spread volatility. We take the SDB stock 
as an example to illustrate the results. The MF-DFA exponent h(q) on q of the 
the adjusted spread return R and the adjusted spread volatility V of the stock 
SDB is shown in Fig. 7 (a) and (d). For a mono-fractal series, one has Ah(q) = 
h(q m in) — h(q max ) = 0. The Ah = h(q_ 6 ) — h(q 6 ) of the spread return and the 
spread volatility is estimated to be 1.02, and 0.15, respectively. It indicates that, the 
spread return shows a strong multifractality, whereas the spread volatility presents 
a weak multifractality. We then investigate the scaling exponent function r(q) on q 
based on partition function, which is also widely adopted to reveal the multifractal- 
ity of a time series, 
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t(q) = qKq) - D f, 



(16) 



where D f is the fractal dimension, with D f = 1 in our case. For the mono-fractal, 
the r(q) linearly increases with q. As shown in Fig. 7 (b), the r(q) of the adjusted 
spread return presents a strong nonlinearity, which further suggests the multifractal 
nature of the spread return. For the adjusted spread volatility, as shown in Fig. 7 (e), 
the r(q) nearly linearly increases with q, which again indicates a weak multifrac- 
tality of the spread volatility. By the Legendre transformation, the local singularity 
exponent ct and its spectrum f(a) can be calculated as pill . 



ct = dr(q)/dq, 
f(a) = qa- r(q), 



(17) 
(18) 




Fig. 7. The MF-DFA exponent h(q) on q is shown in (a) and (d), and the scaling exponent 
function r(q) on q is shown in (b) and (e), and the multifractal singularity spectrum /(a) 
on a is shown in (c) and (f). The upper panels are for the adjusted spread return, and the 
lower panels are for the adjusted spread volatility of the SDB stock, respectively. 



The difference between the maximum and the minimum of the local singularity 
exponent Act = a max — a m i n indicates the width of the extracted multifractal spec- 
trum. The larger the Act, the stronger the multifractality. Fig. 7 (c) and (f) illustrate 
the multifractal singularity spectra f(a) of the adjusted spread return R and the 
adjusted spread volatility V. The width of the extracted multifractal spectrum Act 
is measured to be 1.32 for the spread return R, indicating a strong multifractality, 
and 0.31 for the spread volatility V, indicating a weak multifractality. The results 
of f(a) further support the results obtained from the MF-DFA exponents h(q) and 
the scaling exponent function r(g), i.e., the spread return shows a strong multifrac- 
tality, whereas the spread volatility presents a weak multifractality. 
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6 Conclusion 



In finance, one of the central tasks is to evaluate the market efficiency. The bid-ask 
spread is served as an important indictor to quantify the financial market liquidity 
and vitality. In the present article, based on the tick-by-tick data of four liquid stocks 
in the Chinese stock market, we investigate the memory effect and the multifractal 
nature of the spread return and the spread volatility. It is observed that, the spread 
volatility shows a significant long-range time-correlation, while the spread return is 
short-range time-correlated. Moreover, the spread return is found to exhibit a strong 
multifractal characteristic, which has been observed in dynamics of a number of fi- 
nancial quantities. Differently from the spread return, the spread volatility presents 
a weak multifractal feature. Our work might help to better understand the financial 
market. 
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